o 
o 

(N 
> 

o 

(N 
(N 



Semiclassical Distorted Wave Model Analysis of the (tt , 

Spectrum 



S Formation Inclusive 



M. Kohno,^ Y. Fujiwara,^ Y. Watanabe,^ K. Ogata'* and M. Kawai'' 

^Physics Division, Kyushu Dental College, Kitakyushu 803-8580, Japan 
^Department of Physics, Kyoto University, Kyoto 606-8502, Japan 
'^Department of Advanced Energy Engineering Science, 
Kyushu University, Kasuga, Fukuoka 816-8580, Japan 
"^Department of Physics, Kyushu University, Fukuoka 812-8581, Japan 

{■iv~ , K'^) hyperon production inclusive spectra with = 1.2 GeV/c measured at KEK on ^^C 
and ■^^Si are analyzed by the semiclassical distorted wave model. Single-particle wave functions of 
the target nucleus are treated using Wigner transformation. This method is able to account for the 
energy and angular dependences of the elementary process in nuclear medium without introducing 
the factorization approximation frequently employed. Calculations of the (7r^,if^) A formation 
process, for which there is no free parameter since the A s.p. potential is known, demonstrate that 
the present model is useful to describe inclusive spectra. It is shown that in order to account for the 
experimental data of the E~ formation spectra a repulsive E-nucleus potential is necessary whose 
magnitude is not so strong as around 100 MeV previously suggested. 
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I. INTRODUCTION 

Various meson production reactions in nuclei are a rich 
source of our understanding of hadronic interactions. In 
particular, interactions involving hyperons have to be ex- 
plored by strangeness exchange reactions, since hyperons 
are absent in ordinary nuclear systems. Naturally, the 
absence itself is a consequence of the properties of strange 
hadrons. The interaction between the lambda hyperon 
and the nucleon is fairly well known, since the experi- 
mental data for A hypernuclei has been accumulated in 
last more than 30 years. The T,-N and S-iV interactions, 
by contrast, have not been well understood. Even the 
sign of the S single-particle (s.p.) potential in nuclear 
medium, which reflects basic properties of the Y,-N in- 
teraction, has not been established. Interactions among 
hyperons are far less investigated, except for the AA case. 

In recent years, much experimental effort has been di- 
rected to the study of strange baryons and baryonic res- 
onances in nuclear medium, using incident tt, K and 7 
beams with the energy of 1 '--^ 2 GeV. The extraction 
of meaningful understanding of these baryon properties 
from such reaction processes is not so simple, however. 
For analyzing experimental data we need to take into 
account various effects, such as the proper treatment of 
projectiles and outgoing hadrons, the model description 
of elementary processes in nuclei and the decent descrip- 
tion of the target nucleus and the residual (hyper) nu- 
cleus. We also have to keep in mind the possibility of the 
change of properties of the relevant hadrons themselves 
in nuclear medium. 

Since fully microscopic description is far from practi- 
cal, various approximations are commonly introduced to 
analyze the experimental data. Then, it is important to 
employ a model as simple and reliable as possible, bear- 
ing in mind that lack of some proper treatment can easily 
lead to misunderstanding of the basic hadron properties. 



S formation spectra in (tt, K) and (K, tt) reactions 
with nuclei are not expected to have narrow peaks, be- 
cause of the strong T,N AN coupling. In spite of this, 
however, the early {K, tt) experimental spectra Q were 
interpreted as indicating an attractive S s.p. potential 
with the depth of about 10 MeV The experimental 

discovery of |He [1, has shown that the E-iV interac- 
tion in the T = 1/2 channel is sufficiently attractive to 
support the bound state in this specific nucleus, as dis- 
cussed by Harada It has been recognized, however, 
that due to the strong repulsion in the isospin T = 3/2 
channel, S bound states are unlikely to be observed in 
heavier nuclei. This conjecture was supported by exper- 
imental results on targets of ^Li and ^Be measured at 
BNL f?! . The analysis of the (K, tt) spectra on ^Be from 
BNL [8] given by Dabrowski [9] in a plane wave impulse 
approximation method suggested that the E potential is 
repulsive of the order of 20 MeV. 

The shift and the width of E~ atomic states are an- 
other source of the information on the S-nucleus poten- 
tial. Batty, Friedman and Gal [l^ reexamined the T,~ 
atomic data and concluded that the S potential should 
be attractive at the surface region but changes its sign 
to become repulsive at the higher density region in a nu- 
cleus. 

Theoretical studies for the two-body S-iV force have 
also been inconclusive. In the 1970s, the Nijmegen group 
started to construct hard-core hyperon-nucleon poten- 
tials in a one-boson exchange model. Parameter sets 
corresponding to two typical choices of the SU(3) mixing 
angles were named as models D and F [ll| . The G-matrix 
calculation by Yamamoto and Bando in Ref. [l3| showed 
that the model D yields —29.3 MeV for the E poten- 
tial in nuclear matter at the normal density {kp = 1-35 
fm~^) and the model F repulsive 5.8 MeV, though the 
explicit numbers vary in a different calculational scheme. 
The soft-core versions subsequently constructed by the 
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Nijmegen group [Ts'l tend to predict an attractive S s.p. 
potential in nuclear medium; —27.1 McV in Ref. and 
— 15.3 MeV in the nuclear matter calculation by Schulze 
et al.^. 

A different approach using a non-relativistic SU(6) 
quark model has been developed by the Kyoto-Niigata 
group [13] to obtain a unified description of 

octet baryon-baryon interactions. In this model, the 
description of the short-ranged part of baryon-baryon 
interactions basically provided by the resonating-group 
method with the spin-flavor SU(6) quark model wave 
functions and the one-gluon exchange Fermi-Breit inter- 
action is supplemented by effective meson-exchange po- 
tentials acting between quarks. This model has little am- 
biguities in the hyperon-nucleon sector after the nucleon- 
nucleon interaction is determined. G matrix calculations 
in the lowest order Brueckner theory with the poten- 
tial named as FSS show that the S s.p. potential 
in symmetric nuclear matter is repulsive of the order of 20 
MeV at normal density. The repulsion due to a strongly 
repulsive character in the isospin T = § channel orig- 
inating from quark Pauli effects overcomes an attractive 
contribution in the T = 1/2 ^S'l channel which is similar 
to that in the AA^ case. The latest version of this quark 
model potential, fss2 [13], gives a smaller repulsion of 
about 8 MeV in symmetric nuclear matter. 

We briefly mention relativistic mean filed model de- 
scription for the hyperon sector. This model does not 
seem to have much predictive power, but once param- 
eters are determined to fit basic properties it has wide 
applicability, for example, to a variety of neutron star 
matter calculations. In early models including only a 
meson, the S hyperon is predicted to have the similar 
attractive potential to the A in nuclear medium. Having 
recognized that the E s.p. potential may be repulsive in 
nuclei, the model was extended to include a* meson to 
account for that property. The repulsion of 30 MeV has 
been tentatively used in literature [l^, [13] , although this 
specific number did not have a solid basis. 

It is also noted that Kaiser |21j calculated the S mean- 
field in symmetric nuclear matter in the framework of 
SU(3) chiral perturbation theory and found a moderately 
repulsive potential, that is, 59 MeV for the real part and 
-21.5 MeV for the imaginary one at normal density. 

Recently, inclusive (tt~,K^) spectra corresponding to 
E formation were measured at KEK 22, 23] with better 
accuracy than before, using the pion beam with the mo- 
mentum of = 1.2 GeV/c on medium-to heavy nuclear 
targets. DWIA analyses in Ref. [l^] for ^^Si and similar 
analyses later on other nuclear targets [1^ gave a notable 
conclusion that the E potential is strongly repulsive, as 



large as 100 MeV. Harada and Hirabayashi [2j] showed 
in similar calculations with their optimal Fermi-averaging 
for the elementary t-matrix that the E potential is repul- 
sive inside the nuclear surface, though the actual strength 
varies with the imaginary part supposed. 

The determination of the E-A^ interaction is of fun- 
damental importance in the study of such problems as 
those of neutron star matter and heavy ion collisions, 
because the baryonic component of such hadronic mat- 
ter, especially the hyperon admixture, is governed by the 
basic baryon-baryon interactions. Considering the im- 
portance of determining the Yi-N interaction on the ba- 
sis of experimental data, it is desirable to analyze the 
KEK experiments in a different and independent calcu- 
lational scheme from those in Refs. [12, [13] ■ In this 
paper, we present a semiclassical method for the DWIA 
approach and apply it to (tt^, K^) inclusive spectra. The 
preparatory version of this approach was reported in Ref. 
[25]. The semiclassical distorted wave (SCDW) model 
was originally considered for describing intermediate en- 
ergy nucleon inelastic reactions on nuclei [26| . Applica- 
tions to various {p,p') and (p, n) inclusive spectra |27l.[2^ 
have demonstrated that the method is quantitatively re- 
liable and thus the applications to the wide range of nu- 
clear reactions are promising. 

In Sec. II, we show basic expressions of the semiclassi- 
cal distorted wave model for describing the (tt, K) inclu- 
sive spectra. The formulation using the Wigner transfor- 
mation for the nuclear density matrix is explained. Ac- 
tual optical potential model parameters used for incident 
pions and outgoing kaons are given in Sec. III. Then, nu- 
merical results for the (tt^, AT^) spectra are presented in 
Sec. IV: first for the A formation to see the applicability 
of the SCDW model and next for the S formation. The 
latter case is the main concern of the present paper to 
obtain more solid information about the strength of the 
E single-particle potential than before. In this paper, 
we are concerned with the spectra on light nuclei with 
N = Z, namely ^^Si and ^^C. Conclusions are given in 
Sec. V, with the outlook of the future extension of our 
model. 



II. SEMICLASSICAL DISTORTED WAVE 
MODEL DESCRIPTION OF THE {-k,K) 
INCLUSIVE SPECTRA 

The starting formula for the double differential cross 
section in a standard distorted wave model description of 
the (tt, K) hyperon (Y) production inclusive reaction is 
expressed as 



p,h 

xct^l{r)Mr)M'r'Wh{r')5{W - + th)e{tF - e„), (1) 
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where Xi^'^ ^^'^ X/ ^ represent the incident pion and fi- 
nal kaon wave functions with energies uji and uj f , respec- 
tively, and W — uJi — ujf is the energy transfer. The 
formula describes the process in which the nucleon in 
the occupied single-particle state h is converted to the 
unobserved outgoing hyperon (A or S) state p. The el- 
ementary amplitude of the process + N ^ K + Y \s 
denoted by Vf^p^i^h, which depends on the energy and mo- 
mentum of the particles in the reaction. In order to treat 
such dependence, it is necessary to introduce momentum 
space integration. In that case, the explicit calculations 
involve higher-dimensional integrations. It is desirable, 
in practice, to develop a tractable and trustful approxi- 
mation method. One procedure that has been frequently 
used is the factorization approximation, in which the el- 
ementary process is taken out of the integration, assum- 
ing some averaging wisdom. As used in Ref. [2^ for 
TT~ + p -\- E~, the elementary amplitude in the 

integrand may be replaced by the averaged differential 
cross section over the nucleon momentum distribution 



R. 



da{TT-p K+E-) _ J p{k)%{nk)5{k - P)dk 



J p(k)S{k - P)dk 



(2) 



with P = kx + ky — k^, and is taken outside of the inte- 
gration. The remaining quantity is the Green function, 
which is not difficult to evaluate in the case of a local 
optical potential. A more sophisticated Fermi-averaging 
method was used in Ref. [2J]- Though such procedure 
has been widely applied to show various successes, the 
justification is far from trivial. Important dynamical ef- 
fects might be hidden in the averaging treatment. 



X^Ur±1-) - e±^^-fe-(^)x^^W, (3) 
R±ls) ~ e^''i^-''f'^^hy\R). (4) 



The local momentum k{R) is defined as follows. The di- 
rection is specified by the quantum mechanical momen- 
tum density kq{R) calculated by 



kg{R) 



jR{x<^)*(fl)H)VxW(fi)} 



(5) 



where 3? represents taking the real part, and the mag- 
nitude is determined by the energy-momentum relation 
^k'^{R) + Ur{R) = E at R. Here, Ur{R) is the real 
part of an optical potential for the distorted wave func- 
tion X with energy E. The relativistic energy-momentum 
relation is used for the distorted wave function described 
by the Klein-Gordon equation. 

The above approximation is expected to work well if 
the dominant contributions in the integration over r and 
r' in Eq. (1) is restricted in the region where r and 
r' are close to each other. Actually, the density matrix 
J2h ^*h{'''')4'h{r) brings about this desirable feature, as is 
shown in the following heuristic argument. It is sufficient 
for the qualitative discussion to assume that nuclear s.p. 
wave functions are harmonic oscillator ones. The sum- 
mation over the z-component of the angular momentum 
of each orbit means that we are treating two oscillator 
functions coupled to the total angular momentum L = 0; 



A. SCDW method 

In Ref. [2^, we presented our SCDW approximation 
method for the DWIA cross section formula, Eq. (1). 
There, we introduced a local Fermi gas approximation 
for the target nucleus. In this paper, we improve the 
description by explicitly treating s.p. wave functions of 
the target nucleus. 

The semiclassical treatment was first introduced in the 
description of the intermediate energy nucleon reactions 
on nuclei (26j . Since the amount of numerical calcula- 
tions is reduced, it becomes feasible to include and as- 
sess multi-step contributions. The calculations of (p,_p) 
and (p, n) inclusive spectra have shown that the SCDW 
method works well. 

The semiclassical approximation employs the following 
idea for the propagation of the wave function. Denoting 
the midpoint and the relative coordinates of r and r' in 
Eq. (1) hy R — ^—5^ and s = r' — r, respectively, we 
assume that the propagation of the distorted waves, Xi 
and x/j from i? to r or r' is described by a plane wave 
with the local classical momentum k{R) at the position 



V2£^\nh£h{ri),nh£hir2);L = 0). (6) 

The transformation to the R and s coordinates are car- 
ried out using the Talmi-Moshinsky brackets. 

\nhih{ri),nhih{r2);L = 0) = 

J2 (ni,m;0\nhih,nh£h;0)\neis),NeiR);L ^ 0). (7) 

N,n,e 

The reaction processes which we consider take place 
mostly at the surface of the target nucleus. When R 
is located in the surface region, the dominant compo- 
nents in the right hand side of Eq. (7) are those in which 
2N + £ is the largest. This indicates that the dependence 
on the relative coordinate s is governed by the Os (n = 
and £ — 0) function, which is certainly short-ranged com- 
pared with the size of the target nucleus. Thus we expect 
that the SCDW treatment of Eqs. (3,4) in Eq. (1) is 
meaningful. 

Note that since the SCDW approximation should be 
exact in homogeneous matter, the SCDW works well in- 
side of the nucleus. The above reasoning implies that 
the SCDW approximation is also applicable to the sur- 
face region. This fact is probably connected to the fact 
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that the local density approximation based on the den- 
sity matrix expansion method 29] works well in nuclear 
structure calculations, including the surface region. 



B. Wigner transformation 

In the preparatory calculations in Ref. [2^, we in- 
troduced a Thomas-Fermi approximation for the density 
matrix J^h 4'h{'''')4>h{i') of the target nucleus. Here, we 
elaborate the description of the density matrix by using 
a Wigner transformation. 

The Wigner transformation of the density matrix of 
the target nuclear wave function is defined as 



zK s 



(8) 



^h{R, K) is given by the inverse transformation as 



<^h{R,K) = 



1 



(27r)3 



ds 



xrhiR-ls)MR+ls) 



(9) 



The summation over the z-component of the angular mo- 
mentum is implicit in these expressions. As is shown in 
the Appendix, $/i(i2, may be expressed in terms of 
the Legendre expansion: 

<i>hiR,K)^ J2 Pi{coBKR)<^i{R,K), (10) 
«=even 



where KR denotes the angle between two vectors K and 
R. It is easy to see that the Thomas-Fermi approxima- 
tion for the density matrix used in Ref. [2^ amounts to 
the replacement 



^^P,(cosi?12) ¥^{R,K) 
2 



h e 



-0{kF{R)-K), 



(11) 



where 9{x) is a step function and kpiR) = [3tt^ pr{R)]^^^ 
is a local Fermi momentum determined by the local pro- 
ton or neutron density Pt{R) at R. 

We describe unobserved hyperons by a local optical 
potential. Actually, the hyperon optical potential should 
be complex, because there are inelastic processes. The 
standard way to treat the completeness of the final states 
described by a complex Hamiltonian is to use the Green 
function method 30]. At this stage, however, we adopt a 
simplified prescription, employing a real potential of the 
standard Woods-Saxon form, 

C/y(r) -;7^/{l + exp((r-ro)/a)}, (12) 
and we convolute the result of the calculated spectrum 
with a Lorentz-type distribution function with the typical 
width, simulating the effects of inelastic channels. In that 
case, the expression of Eq. (1) is directly used. Then, we 
introduce the SCDW approximation as in Eqs. (3,4) also 
for the hyperon wave functions 4>p{r) and (f>p{r'). 

Using these approximations explained above, the dou- 
ble differential cross section in the SCDW model becomes 



a 



fdRfdKY: --^lx^-^(H)l^lx!+)(i2)1^10p(H)l^(27r)3^^P,(cosKi?)$f,(i?,i^) 



dWdn (27 

X \vf,p,,,h{K, k,)\-'S{K + k,{R) - kf{R) - kp{R)) 6{W - + e,,) 

I 



(13) 



The summation means the sum over the spin and 
the momentum of the outgoing unobserved hyperon: 
■(2^ X^spin / '^P- ^-'^^ energy, instead of the 

momentum, to specify scattering states, the momentum 
integration J dp is written as follows: 

/ ^^p/K^)'^'^"^'^- ^^^^ 

The above final expression (|13p admits of the simple 
interpretation that the reaction in which tt+N yields K+ 
Y takes place at the position R and satisfies conservation 
of local semiclassical momentum: 

K + k,{R) ^ kf{R) + kp{R) (15) 



These momenta, ki (R) , k f {R) and kp {R) , are calculated 
with Eq. (5), using tt, K and Y distorted wave functions 
in an optical model description. It is seen that the di- 
mension of the integration does not change from Eq. (1) 
to Eq. (jl3p . However, we can now treat the momentum 
dependence of the transition amplitude w/,p,i,;i explicitly. 

In the present formulation, the correction for the lack 
of the translational invariance [31I ] in the target nuclear 
wave function is not included. In describing electron 
scattering on a nucleus composed of A nucleons, the 
center-of-mass effects 32] in the nuclear shell model have 
been taken care of by a multiplicative factor F^/^ = 
ex.p{q'^ / (AAa)) , where q is the momentum transfer and 
a is the oscillator constant. For A = 12 and a — 0.4 
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fm~^, the factor amounts to 1.5 when the mo- 

mentum transfer is 2 fm~^. Thus the calculation without 
the center-of-mass correction tends to underestimate the 
cross section. The appropriate treatment of the center-of- 
mass effects in our SCDW method deserves to be studied 
in the future. 



C. Elementary amplitude 

The on-shell amplitude Vf^p^i^h of the elementary pro- 
cess is related to the differential cross section by 



da 



1 



EnEy kx . 1 2 
\v\ , 



(16) 



where s is the invariant mass squared. In Eq. p^ . we 
are able to account for the angular dependence of the 
TT + N K + Y elementary process at each position 
R. To carry out an actual calculation in Eq. ^3]) . we 
need some model description for v. However, at present, 
there is no reliable model for the relevant process, includ- 
ing off-shell regions. We use a simple phenomenological 
parametrization based on Eq. (16), by defining the in- 
variant mass squared using the momenta ki{R) and K. 
The following functional form is used to simulate empiri- 
cal angular distributions of the tt + N K + Y reactions. 



da 

dn 



An 



{l + J2ai{V^)Pe{coae)} (17) 



Values of a{y/s) and ai(yG) used in the present calcula- 
tions will be given in Sec. IV. 



D. Wave functions of the target nucleus and 
hyperons 

In this paper, we are concerned with the reactions on 
the N = Z targets, ^*Si and ^^C. Single-particle wave 
functions and the energies for these nuclei eu are pre- 
pared by the density-dependent Hartree-Fock description 
of Campi-Sprung |33| . 

As stated in B of this section, hyperons are described 
by an energy-independent local potential of the Woods- 
Saxon form. We use the standard geometry parameters, 
ro = 1.25 X {A- 1)1/3 fm and a = 0.65 fm. The Coulomb 
potential regularized in a nucleus is incorporated in the 
case of the S^. It is noted that if we use a different 
parameter set such as rg = 1.20 x {A — 1)^/'^ fm and a — 
0.6 fm, we do not see appreciable changes in calculated 
spectra for ^^C and ^^Si. 

It may be argued that the hyperon s.p. potential is not 
expressed by a single Woods-Saxon shape. The S poten- 
tial may be repulsive in the higher density region, but 
change its sign at the surface as has been suggested by the 
analysis [l3] of S~ atomic data and also by the nuclear 
matter calculations with the SU(6) quark model in- 
teraction. At present, however, it is premature to discuss 
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FIG. 1: Differential cross section of tt elastic scattering on 
^^C at Ptt = 800 MeV/c. Calculated values in the optical 
potential model with — i6o ~ 1-0 fm^ is compared with the 
experimental data [34| . 



the detailed shape of the E potential from the available 
inclusive spectra. We assume, from the beginning, the 
standard Woods-Saxon form and ask what strength C/£ 
is favored by the experimental data. More experimen- 
tal data is needed to quantitatively discuss more elab- 
orate shape parameters as well as the possible energy- 
dependence of the strength. 

The nuclear matter calculations [l^ with the SU(6) 
quark model interaction suggest that the imaginary 
strength of the A s.p. potential hardly depends on the 
density, but that of the S s.p. potential decreases when 
the nuclear matter density becomes lower. Regarding 
that the hyperon formation processes take places at the 
surface region, we set the energy dependence of the half 
width in MeV for smearing the calculated spectra as fol- 
lows, simulating the calculated results at the half of the 
normal density with the quark model potential ESS: 



^a{E) 
2 



1 + 12(^/40)2 -8(£;/40)^ S<40MeV 

(18) 

5, £; > 40 MeV, 
= 10+^ £;inMeV. (19) 



1 + (20/£;)2' 



III. OPTICAL POTENTIALS FOR PIONS AND 
KAONS 

The incident pions and detected kaons are described 
by the standard Klein-Gordon equation with some op- 
tical potential model. Following the usual procedure to 
construct a 7r-nucleus optical potential from ttN elemen- 
tary amplitudes, the optical potential for the pion is given 

by 



(20) 
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FIG. 2; Differential cross section of elastic scattering on 
^^C at pk = 800 MeV/c. The optical model with the param- 
eter 6o given in Eq. (21) gives the result shown by a solid 
curve, compared with the experimental data [3^ . 




cm. angle [deg] 

FIG. 3; Differential cross section of elastic scattering on 
^^C at Pk ~ 715 MeV/c. The optical model with the param- 
eter ho given in Eq. (21) gives the result shown by a solid 
curve, compared with the experimental data [37| . 



TABLE I: Legendre coefficients ai{^) in Eq. (17) of the 
differential cross sections of the tt" -\-p K'^ + A reaction, 
determined using experimental data [iol . l4l| . 





ai range of -^s (GeV) 


e = i 


12.6Vs - 20.362 < 1.69 

0.932 1.69 <^< 1.925 
— 4. 8v^-l- 10.172 1.925 < ^ < 2.015 
2.164^5- 3.86 2.015 < < 2.3 


1 = 2 


Vs< 1.677 
11.82^5- 19.821 1.677 <ys< 1.7 
0.273 1.7 < ^ < 1.8 
4.71^3- 8.205 1.8 <^<2.4 


£ = 3 


Vs<1.8 
-5.68Vs + 10.224 1.8 <Vs< 1-85 
-0.284 1.85 <ys< 1.98 
-8.32^5 + 16.19 1.98 <Vs< 2.03 
6.13^5- 13.144 2.03 <^< 2.4 


£ = 4 


^/s<1.86 
7.5^7- 13.95 1.86 <^< 1.96 
0.75 1.96 < ^/s < 2.08 
5.5^/s- 10.69 2.08 < < 2.4 



search an optimal parameter 5o in the form of Eq. (20), 
using the available experimental data [sl, H^l on C at 
800 MeV/c and 715 MeV/c. For the K+, we find that 
the following momentum dependence is adequate; 

on 

6o = -6.2 X 10" V' + i— , (21) 

VK 

where 6o in fni"^ and "pK in MeV/c. As shown in Figs. 2 
and 3, the calculated differential cross sections account 
well for the experimental data. The outgoing K'^ mo- 
menta relevant to the present inclusive spectra arc in this 
energy region. 

IV. RESULTS 
A. A formation 



where p{r) is the one-body nuclear density distribution 
and the parameter 6o is related to the sum of isospin av- 
eraged ttTV partial wave amplitudes. In practice, a pure 
imaginary choice of &o = ^\{'^tot) is found to work well. 
As an example, we show, in Fig. 1 , the pion elastic scat- 
tering differential cross sections on ^^C at 800 MeV/c. 
We simply expect at the present stage that the same 
prescription is applicable to the incident momentum of 
1.2 GeV/c. Actually, -ih^ = 0.58 fm^ at 1.2 GeV/c and 
-ihQ = 0.70 fm3 at 1.05 GeV/c. 

It has been known _35] that the kaon scattering data is 
not well reproduced by simply folding elementary ampli- 
tudes. In the present calculation, we phenomenologically 



We first apply our SCDW model to the [11+, K+) A 
formation inclusive spectra obtained with ^^Si and ^^C 
targets measured at KEK [H,!!!]. Since the A s.p. po- 
tential has been established as ~ —30 MeV from var- 
ious A hypernuclear data ^3^, the calculation serves as 
quantitative assessment of the SCDW model description. 
The strength and angular dependence of the elementary 
process are parameterized as Eq. (17) accordin g to the 
available experimental data of tt^ + p ^ + A |40l. l4lj|. 
The energy dependence of the total cross section is fitted 
as 



1.2 



472(yi- 1.60935) 
[26(Vs- 1.60935)]3 20 



[mb]. 



(22) 
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FIG. 4: (tt"*", Ti'"'") A formation inclusive spectra with a ^*Si 
target at 9k ~ G° T 2° for pions with = 1.2 GeV/c. These 
results were obtained with three choices of U'^ in a Woods- 
Saxon potential form with the geometry parameters of ro = 
1.25 x{A- 1)^/3 fm and a = 0.65 fm. The KEK data [H are 
also displayed. 



FIG. 5: (tt"*", is'"'") A formation inclusive spectra with a ^'^C 
target at 6k = 6° T 2° for pions with p-^ = 1.2 GeV/c. These 
results were obtained with three choices of U'^ in a Woods- 
Saxon potential form with the geometry parameters of ro = 
1.25 X (A - 1)^/^ fm and a = 0.65 fm. The KEK data 
are also displayed. 



Coefficients 0^(^/5) for the angular dependence are tab- 
ulated in Table I. 

Calculated spectra with three choices of the strength of 
the A s.p. potential, = -40, -30, -10 MeV, are com- 
pared to see the potential dependence. Figure 4 shows 
the results for ^*Si with p^^ — 1.2 GeV/c, while Figs. 5 
and 6 for ^^C with = 1.2 and 1.05 GeV/c, respectively. 

The A hyperon can be bound in a nucleus. In ^^C, Os 
and Op A bound states appear at energies of —12.4 MeV 
and —1.8 MeV, respectively, for the case of = —30 
MeV with tq = 1.25 x ll^/^ fm and a = 0.65 fm, ne- 
glecting the small spin-orbit component. The A bound 
state wave function is treated in the same manner as for 
the target nuclear wave function, Eq. (9). The scat- 
tering wave function (/)p in Eq. (13) is replaced by the 
Wigner transformation of the A bound state wave func- 
tion. The transition strength from the 0^3/2 nucleon hole 
state appears as two peaks below — Sa = 0, for which 
the experimental resolution of 2 MeV in FWHM is ap- 
plied. In this presentation, we use s.p. neutron energies 
in a simple 0p3/2-closed HE description for the target 
nucleus ^^C. Thus the peak position does not precisely 
agree with that of the experimental spectrum. Since the 
Osi/2 nucleon hole state has a large width, we broaden 
the calculated strength by the half width of 8 MeV. 

These figures indicate that the energy dependence of 
the inclusive spectrum is well reproduced by the standard 
choice of the A s.p. potential, Vq = —30 MeV. The abso- 
lute values are short by about 35 %. As is noted in the 
end of Sec. II B, our calculation tends to underestimate 
the cross section for lack of the translational invariance in 
the target wave function. We also expect various other 
effects for the underestimation. Besides possible ambi- 
guities in the SCDW treatment and uncertainties in the 



elementary strengths as well as the optical potential pa- 
rameters, there should be room for contributions from 
two-step processes and more. Since the incident pion ab- 
sorption is rather large, some of the flux lost may emerge 
again into the A production channel. Bearing in mind 
these points and in addition the possible modification 
of the elementary process in nuclear medium, which is 
clearly premature to discuss at the present stage of the 
analysis, our SCDW model is considered to provide a 
meaningful description for the (tt, K) inclusive spectra. 
The explicit estimation of the center of mass correction 
and the two-step contributions is needed to establish the 
quantitative reliability of the SCDW method. 



B. E formation 



In Ref. 25|, we assumed an isotropic angular depen- 
dence for the 7r~ + p — > + T,~ elementary process 
in calculating (7r~,i^+) formation inclusive spectra. 
The energy dependence of o'{^/s) was taken from the pa- 
rameterization by Tsushima et al. 0| , which was renor- 
malized by a factor of 0.82. In this paper, we take into 
account the angular distribution using the Legendrepoly- 
nomial coefficients reported by Good and Kofler [43| on 
the basis of available data [H, H, Hi, Hi] . Up to ~ 2.1 
GeV, - 1900 MeV/c, we can set = with £ > 3. 
We try to simulate the energy dependence of ai(\/s) and 
02 (Vs) by several lines as given in Table II, which are de- 
picted in Fig. 7 with the experimental data. The energy 
dependence of the total cross section <j{\fs) is parame- 
terized as follows: 



0.9603 



0.02055(Vi- 1.691) 
(Vs- 1.682)2 -^0.003131 
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= 1 .05 GeV/c 




50 100 

-Ba [MeV] 

FIG. 6: {■K'^yK'^) A formation inclusive spectra with a ^^C 
target at 9k = 6° T 2° for pions with p,r ~ 1.05 GeV/c. 
These results were obtained with three choices of in a 
Woods-Saxon potential form with the geometry parameters 
of ro = 1.25 X {A - 1)^/^ fm and a = 0.65 fm. The KEK data 
[SSf l are also displayed. 
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FIG. 7: y^-dependence of Legendre coefficients of angular 
distributions of the ■k~ + p ^ + E~ reaction. Empirical 
data points are taken from Refs. [H, IH, lil, [ii] . The dotted 
and dashed lines show the fit for the calculation of (Tr",^^^) 
inclusive spectra. 
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TABLE II: Legendre coefficients af(y's) in Eq. (17) of the 
differential cross sections of the 7r~ + p + E~ reaction, 

determined using experimental data [431 . |44| . |45| . [i^ . 





ai 


range of ^/s (GeV) 


i = 1 




0.0 


^/s < 1.72 




-7.134(ys- 


- 1.72) 


1.72 <y/s< 1.93 






-1.5 


1.93 < ^< 1.97 




6.25^5- 


- 13.81 


1.97 < ^/s < 2.05 






-1.0 


2.05 <y/s< 2.15 




-4.06Vs H 


- 7.729 


2.15 < < 2.31 






-1.65 


2.31 < y/s 


i = 2 




0.55 


^/i < 1.81 




S.O^s- 


- 13.93 


1.81 < < 1.89 




-7.44^5+15.25 


1.89 <y/s< 2.05 






0.0 


2.05 < ^/s < 2.15 




lo.oys 


- 21.5 


2.15 < < 2.31 






1.6 


2.31 < ^ 




Ql—L , , \ , , 1 

1.7 1.9 2.1 2.3 

/s [GeV] 



FIG. 8: Total 7r~ +p + S~ cross sections as a function 

of the cm. energy. Data points are taken from Refs. [isl . 
El, l45l . H^. The solid curve displays the parameterization 
of Eq. (23) in the present calculation. For comparison, the 
energy dependence used in Ref. with an isotropic angular 
dependence is shown by the dotted curve. 



0.003023(V^- 1.691)°"9^ 
(Vs- 1.894)2 + 0.01548 



[mb], (23) 



which is displayed by the solid curve in Fig. 8. At 
^/s = 1.79 GeV, corresponding to = 1.2 GeV/c, the 
absolute magnitude of the differential cross section at 
forward angles in the laboratory frame is about 1 30 u b, 
which corresponds well to that measured at KEK |23 |. 

Figures 9 and 10 compare calculated (tt", Jtr+) inclu- 
sive spectra with the KEK experimental data [23j . 
Several curves in these figures correspond to the assumed 
S potential with (75 = -10, 10, 30 and 50 MeV. In DWIA 
analyses in Refs. [l^, and also in Ref. [13] , an overall 
reduction factor is introduced to discuss the correspon- 



dence with the experimental data. However, we do not 
multiply any renormalization factor. It is seen that ab- 
solute values are satisfactorily reproduced by a repulsive 
strength of 10 ~ 30 MeV. Since we may expect addi- 
tional contributions from multi-step processes, the actual 
repulsive strength may be larger than 30 MeV. This re- 
sult agrees with that in Ref. [2g| , where the S potential 
was concluded to be repulsive of the order of 30 ~ 50 
MeV, using the SCDW model with the Thomas-Fermi 
approximation for the target nucleus, ^^Si. The assump- 
tion of the isotropic angular distribution of the elemen- 
tary process used in Ref. [25] tends to overestimate the 
elementary cross section at forward angles. Thus, the 
repulsive strength needed to reproduce the experimental 
data becomes smaller. 



9 



> 

CD 



t 2 



J3 
LU 

a 

"a 

T3 



'°Si(7i",K") 
p^- = 1 .2 GeV/c 



Ur 

-1 MeV- 
+1 MeV. 
+30 MeV- 
+50 MeV- 




100 200 
-Bj; [MeV] 



FIG. 9: (Tr",^^ ) S formation inclusive spectra with a Si 
target at 9k = 6° =F 2° for pions with p-^ = 1.2 GeV/c. 
These results were obtained with four choices of the repul- 
sive strength = -10, 10, 30, 50 in a Woods-Saxon potential 
form with the geometry parameters of ro = 1.25 x {A — 1)^''^ 
fm and a=0.65 fm. Experimental data points are taken from 

Refs. in,!!!. 
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FIG. 10: {tv' ,K'^) E formation inclusive spectra with a ^^C 
target at 9k = 6° =F 2° for pions with p,r = 1.2 GeV/c. 
These results were obtained with four choices of the repul- 
sive strength = -10, 10, 30, 50 in a Woods-S axon potential 
form with the geometry parameters of ro = 1.25 x (A — 1)^'^^ 
fm and a=0.65 fm. Experimental data points are taken from 

Refs. ;23;i. 



At present, the agreement of the calculated shape of 
the spectrum with the experimental data is not excel- 
lent. This may be related to the uncertainties of the in- 
put cross section, besides multi-step contributions. As is 
seen in Figs. 7 and 8, error bars of the elementary cross 
section are rather large. It is also to be born in mind 
that the S-nucleus potential may be energy-dependent. 
On the experimental side, the data is uncertain at the 
higher excitation energy region due to the spectrometer 



TABLE TIL Percentage contribution from each cm. energy 
(v^) region of the ■n~ -\-p ^ + E~ elementary process at 
-Be = 50, 110 and 170 MeV in {tt',K+) E formation inclu- 
sive spectra on ^®Si with — 10 and 30 MeV, respectively. 



case of 


= 10 MeV 






-Bs (MeV) 


range of a/s (GeV) 


50 110 


170 


1.75 


0.0 % 0.0 % 


3.7 % 


1.75 <^/s< 1.80 


0.0 % 25.7 % 


68.2 % 


1.80 <^/s< 1.85 


21.8 % 62.7 % 


27.3 % 


1.85 <^/s< 1.90 


68.4 % 11.2 % 


0.8 % 


1.90 < ^< 2.00 


9.7 % 0.4 % 


0.0 % 


2.00 < ^/s 


0.1 % 0.0 % 


0.0 % 


case of 


= 30 MeV 






-Be (MeV) 


range of -y/s (GeV) 


50 110 


170 


^< 1.75 


0.0 % 0.0 % 


1.8 % 


1.75 <y/s< 1.80 


0.0 % 18.1 % 


48.3 % 


1.80 <^/s< 1.85 


12.2 % 55.2 % 


46.2 % 


1.85 < ^< 1.90 


66.6 % 25.0 % 


3.6 % 


1.90 <y/s< 2.00 


20.9 % 1.7 % 


0.1 % 


2.00 < 


0.3 % 0.0 % 


0.0 % 



acceptance [23|. 

In order to learn the role of the nucleon Fermi mo- 
tion, it is instructive to show which energy region of the 
elementary process dominantly contributes to the forma- 
tion strength at each S separation energy, —B^. Table 
III tabulates percentage contributions from six different 
regions of ^/s at -B^ = 50, 110 and 170 MeV for the 
cases of = 10 and 30 MeV, respectively. At lower 
E excitation energies, the reaction mainly occurs with a 
nucleon moving toward the incident pion. On the other 
hand, at higher S energies, the dominant contribution 
comes from the elementary process at lower cm. ener- 
gies. 

It was remarked in Ref. [l^l that the peak position 
of the {t:~ , K~^) E~ formation inclusive spectrum at an 
energy as high as 120 MeV is difficult to reproduce if the 
repulsion of the E-nucleus potential is not so strong as 
about 100 MeV. Our analysis suggests, however, that it 
is not necessary for the E s.p. potential to be such repul- 
sive. The reason that the result obtained in our SCDW 
model differs from that of Ref. [l^l might be related to 
the fact that we did not use the factorization approxima- 
tion represented by the average cross section, Eq. (2). 
The descriptions of incident pion and outgoing kaon dis- 
torted waves are also different. It will be worthwhile to 
detect the source of the different results. 

It is necessary to discuss the correspondence of the E 
potential strength obtained on the basis of the present 
(Tr~ , K~^) inclusive spectra with the predictions of vary- 
ing theoretical models for the EA^ interaction. Most 
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models [ll|, [iSj of the Nijmegen group give an attrac- 
tive E s.p. potential, as is seen in various G-matrix 
calculations [H, 0, [13, SH] in nuclear matter. Only 
the model F is acceptable, which was concluded also by 
Dabrowski Q in his plane wave impulse approximation 
analysis of the BNL data \§\ of the {K~ , 7r+) spectrum on 
^Be at p k = 600 MeV/c. The recent SU(6) quark model 
[isl . [TgL [i3] by Kyoto-Niigata group definitely predicts 
a repulsive S s.p. potential. The G-matrix calculations 
fisl in symmetric nuclear matter showed that the early 
H, [li, gives t/g - -t-20 MeV and the lat- 
13, J7° ~ -f8 MeV. This repulsive E 



version, FSS 
est version, fss2 



s.p. potential originates from a strongly repulsive char- 
acter in the isospin T = | channel due to the quark 
anti-symmetrization effect. If the strength of more than 
30 MeV is confirmed in future, these theoretical models 
will need fine tuning. 

It is also important to pay attention to the density de- 
pendence of the E s.p. potential. As the E~ atomic data 
tells , the E potential is attractive at the very surface 
region of a nucleus. This feature is also seen in the G- 
matrix calculation with the quark model potential [l8l |. 
The calculation in this paper assumes a single Woods- 
Saxon form for the E potential. A question whether the 
sign change of the E^ potential can be detected or not 
at the outside region in the description of the (tt", if +) 
reaction is to be studied in future. The energy depen- 
dence of the E potential is another issue to be addressed. 
The G-matrix calculation in Ref. ^18] indicates that the 
repulsive strength is not energy-independent. 

Because of the strong repulsive contribution from the 
isospin T = I channel, it is hypothesized that the E~- 
nucleus potential becomes more repulsive in the case of a 
neutron excess. In this respect, analyses of the (tt^, K^) 
data with heavier nucleus targets will be interesting for 
the purpose of investigating whether such quantitative 
isospin dependence actually exists. 

In the present calculations, there are various treat- 
ments to be improved. The smearing caused by the 
Lorentz-type convolution should be treated by the pre- 
cise way of the Green's function method, though much 
calculational efforts have to be devoted. A quantitative 
estimation of the contribution from multi-step processes 
is needed, which could fill the difference of the experimen- 
tal data and the calculational results as is seen in the A 
formation spectra. A model description of the elemen- 
tary process is to be upgraded, though new experimental 
data is required to do it. After improving these points 
and the proper account of the CM motion of the target 
wave function, the SCDW framework serves as a quan- 
titatively reliable model to study the possible modifica- 
tion of the elementary amplitudes in nuclear medium. 
The direction of the change of cross sections, increase 
or decrease, depends on how the amplitudes are altered 
through the underlying dynamical processes. 



V. CONCLUSIONS 

We have developed a semiclassical distorted wave 
model for (tt, K) inclusive spectra corresponding to A 
and E formation processes, using the Wigner transfor- 
mation of the nuclear density matrix. The expression 
of the double differential cross section consists of the in- 
coming pion distorted wave function, the outgoing kaon 
distorted wave function and the undetected hyperon dis- 
torted wave function at each collision point, where the 
conservation of the classical local momenta is respected. 
The momentum distribution of the bound nucleon in the 
target nucleus is obtained from the Wigner transforma- 
tion of Hartree-Fock wave functions. 

We have first applied the model to inclusive (tt"*", K^) 
A formation spectra on the ^*Si and ^^C targets mea- 
sured at KEK [2^. In this case, since the A s.p. po- 
tential is well known, there is no adjustable parameter. 
The standard A potential strength is found to reproduce 
well overall energy dependence of the data. The strength 
is underestimated. However, this is a rather preferable 
result, because the proper treatment of the CM motion 
has to be implemented in our SCDW formulation and 
there should be some contributions from two-step pro- 
cesses which are not taken into account. The quantita- 
tive estimation of these effects is one of the important 
subjects to be investigated. 

Observing from the (7r+, if +) A formation spectra that 
the SCDW model provides a useful description of the 
inclusive spectrum without any adjustable parameters 
and renormalization factors, we have proceeded to the 
(Tr~,K~^) E formation spectra. The comparison of the 
calculated curves using several choices of the E s.p. po- 
tential strength in a standard Woods-Saxon geometry 
with experimental data from KEK [23J has shown 
that an attractive E-nucleus potential overestimates the 
spectrum at lower excitation energies. Although there 
are rather large uncertainties in the information of the 
elementary process, we see that the repulsive potential 
is necessary to account for the absolute strength of the 
spectrum. Although we have to await quantitative es- 
timation of various effects above mentioned to specify 
the strength of the E-nucleus potential, it is reasonable 
to conclude that the E hyperon experiences repulsion in 
nuclear medium and its magnitude is not so strong as 
around 100 MeV which was suggested by DWIA analy- 
ses in Ref. [13 • 

The information about the repulsive feature of the E- 
nucleus potential constrains the two-body E-A^ poten- 
tial model and thereby improves our understanding of 
the interactions between octet baryons. In the literature 
there has been a few E-A^ potential models which 
diet repulsive E mean field. In Nijmegen models [ll|, 
only the model F is satisfactory in this respect. An- 
other model is a SU(6) quark model by the Kyoto-Niigata 
group 0,[li,[l3. The model FSS gives 20 MeV [11], on 
the other hand the more sophisticated version fss2 pre- 
dicts somewhat smaller repulsion of 8 MeV. More studies 
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are certainly needed to determine the strength of the S 
s.p. potential, by employing various choice of the po- 
tential shapes. The energy dependence of the E-nucleus 
potential may also have to be taken into consideration. 

The analyses of the data of heavier target nuclei are im- 
portant in the next step. Since the neutron excess means 
that the T = | contribution becomes larger, we could 
check the isospin dependence of the E-7V interaction on 
the basis of experimental data. The SCDW analyses of 
the data on ^^Ni, "^In and ^°^Bi taken at KEK j23|] are 
in progress. 

Finally we note that the present framework can be 
straightforwardly extended to describe other inclusive 
spectra, such as {K,tt), {K^,K^), (71,77), {l^K), {j,ri) 
and so on. 

This study is supported by Grants-in-Aid for Scientific 
Research (C) from the Japan Society for the Promotion of 
Science (Grant Nos. 15540284, 15540270 and 17540263). 



APPENDIX A: WIGNER TRANSFORMATION 
OF THE DENSITY MATRIX 

We present an explicit expression for the Wigner trans- 
formation of the density matrix: 



MR, K) 



(27r)3 

xJ2 'I^UR^ls)MR + Is). (Al) 



We write the s.p. wave function of each partial wave in 
r-space as 

Mr) = lM,e,.j,Ar)[yi,{r) x XiA, (A2) 

where (t'nh,ih,jhi''') ^ radial wave function and X1/2 ^ 
spin part. Let us denote the Fourier transform of the 
single-particle wave function fphir) as 4ih{k). 

Mk) = / dre-'^-^'Mr) 



(27r)3 



{2Ti)l 
X-^^nhAjh(fc)> 



(A3) 



where k represents angular parts of the vector k and 
the Fourier transformation of the radial wave function is 
defined as 

^<^nhAjh(*) = (-*)^"''v/|- / drrj/,^{kr)(l}„^j^^j^{r). 



^hiR,K) in Eq. (Al) is rewritten as 



(A4) 



MR,K)^ I dpMK-\p)MK+]:P)^P-^. (A5) 



Using the expression of Eq. (A3), we first obtain 
Y,MK ^\p)k{K + \p) 

rrih 

AYiAK + ]:P) y. X1/2] 



\K~\p\ \K+\p\ 



-\p\)M.U.,u{\K +\p\). (A6) 



Since the spin part gives m'^ ^ ^m^m'^, the 

recoupling of angular momenta leads to the following ex- 
pression. 



Y,4>liK-lp)kiK + lp) 



1 



= (^^^^^(^°^'^^P) \K - ipl \K + ipl 

ip|)^„„,,.,,(|X + ip|). (A7) 

Here, the angle between K — ^p and K + ^p is denoted 
by p] that is, 



cos^ 



(A8) 



Then, Eq. (A5) becomes 

MR^K). jdp,np^^J^ 

X p^^^^^^^k^p)\k-Ip\\kIip\ 

X ~^1^^,^^^^{\K -\p\)i^,,,,,,{\K +\p\). 

(A9) 

Noting that the following relation holds 

r.27r 







Rp 



27:J2^\2e+l)MRp) 



X PeicosKR)Pe{cosKp),{A10) 



we obtain the Legendre expansion of ^ii{R, K). Each 
component $^(_R, AT) of the expansion 



MR,K)^ Pi{cosKR)¥y^{R,K) (All) 

is given as follows 



«=even 
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^i{R,K) = 27r(-l)^/2(2£+l) J p^dpdcosKpje{Rp)Pi{cosKp){2jh + l)j^^PedooseKp) 

2 f 1 f 1 

^-Jdr rje^{\K + -p|r)</)„^,^^j^(r) J dr rjt^{\K - ^PlO'/'n^ A W- (A12) 
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